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ABSTRACT 


Isenthalpic  solidification  of  a  pure  supercooled  liquid  Is  shown  to 
result  In  either  a  two-phase  solld/llquid  mixture  In  Invariant  equilib¬ 
rium  or  a  single-phase,  totally  solid  material  In  unlvarlant  equilibrium, 
depending  on  the  level  of  supercooling  prior  to  solidification.  The  crit¬ 
ical  supercooling  above  which  unlvarlant  equilibrium  Is  obtained  is 
large  for  metals  (hundreds  of  centigrade  degrees)  but  much  smaller  for 
certain  molecular  substances.  Experiments  on  white  phosphorus  (al¬ 
pha  P4)  show  that  the  critical  supercooling  (25.6  C°)  can  be  reached, 
and  exceeded,  easily.  Solidification  rate  measurements  taken  above 
and  below  the  critical  supercooling  for  P4  show  that  the  solld/llquid 
interface  temperature  varies  smoothly  with  melt  supercooling,  although 
light-scattering  experiments  Indicate  that  rapid  changes  occur  In  the 
extent  of  the  dendritic  zone  as  the  critical  supercooling  Is  approached 
and  exceeded. 

A  method  for  extracting  interface  attachment  kinetics  from  solidi¬ 
fication  rate  data  was  examined  In  detail  and  applied  to  our  rate  meas¬ 
urements  on  P4.  We  find  that  above  about  9  C3  supercooling,  P,  solidi¬ 
fies  with  linear  attachment  kinetics  having  a  rate  constant  of  17.7  ±  0.4 
cm/sec  C*.  Below  about  1  C°  supercooling,  P4  solidifies  with  a  faceted 
morphology  Indicative  of  layer-passage  limited  kinetics.  Between  1  C° 
and  9  C*  supercooling,  transitional  growth  kinetics  occur.  Thexe  re¬ 
sults  are  in  qualitative  agreement  with  the  crystal  growth  theory  of 
Cahn,  et  al.,  which  predicts  that  attachmert  kinetics  should  change  as 
the  driving  force  for  crystal  growth  is  varied  by  substantial  amounts. 
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INVESTIGATION  OF  SOLID/LIQUID  INTERFACE 
TEMPERATURES  VIA  ISENTHALPIC  SOLIDIFICATION 


INTRODUCTION 

Recently,  several  successful  experimental  approaches  have  been  described  in  which 
slowly  advancing,  complex  dendritic  interfaces  in  freezing  organic  compounds  (1,2)  and 
metals  (3,4)  were  analyzed  on  the  basis  of  in  situ  microscopic  observation.  Because  the 
solidification  speeds  reported  in  those  studies  were  near  the  lower  limit  for  dendritic 
growth  (e.g.,  ~10  2  cm/sec  in  pure  rnetals),  the  results  were  analyzable  solely  on  the 
basis  of  heat-flow  theory  for  isolated  dendrites;  no  consideration  was  given  to  atomic 
attachment  processes  at  the  solid/liquid  interface.  Hamilton  and  Seideusticker  (5)  have 
reported  on  an  excellent  technique  for  estimating  the  growth  kinetics  of  germanium  den¬ 
drites  growing  in  the  speed  range  0.1  to  1  cm/sec.  Their  technique  is  applicable  to  ma¬ 
terials  in  which  single  dendrites  or  ribbon  crystals  can  be  grown  under  near  steady-state 
conditions.  These  successes  notwithstanding,  no  suitable  method  has  yet  been  demon¬ 
strated  for  estimating  the  solidification  kinetics  of  substances  that  solidify  by  rapid  den¬ 
dritic  growth  in  the  speed  range  10  to  10  4  cm/sec.  The  lack  of  success  in  this  area  is 
perhaps  not  surprising  considering  our  ignorance  about  the  morphological  details  of  the 
solid/liquid  interface,  which  is  topographically  complex  and  advancing  at  high  speed. 
Without  this  precise  morphological  information,  a  reliable  estimate  cannot  ordinarily  be 
made  concerning  the  departure  of  interface  temperature  from  the  local  equilibrium  tem¬ 
perature;  consequently,  a  kinetic  relationship  cannot  be  established  between  the  magni¬ 
tude  of  this  temperature  departure,  which  is  the  driving  force  for  the  molecular  proc¬ 
esses  of  solidification,  and  the  speed  of  advance  of  the  interface  normal  to  Itself,  which 
is  the  local  kinetic  rate  of  solidification. 

The  purpose  of  this  report  Is  to  present  an  approach  to  the  problem  of  obtaining  In¬ 
formation  on  the  kinetics  of  rapid  solidification  processes.  Specifically,  an  analysis  of 
solidification  at  extreme  levels  of  supercooling  will  be  presented  which  yields  a  method 
for  establishing  limits  on  the  solid/liquid  interface  temperature  at  the  kinetlcally  most 
active  regions  of  the  interface.  Several  experiments  on  supercooled  white  phosphorus 
will  then  be  described  which  provide  a  basis  for  the  theoretical  analysis  and  which  yield 
quantitative  kinetic  data  on  molecular  attachment  processes  during  rapid  solidification. 


THEORY 
Thermodynamic  s 

Consider  a  uniformly  supercooled  pure  melt  thermally  isolated  from  its  surround¬ 
ings.  The  initial  (metastable)  equilibrium  state  of  this  system  is  characterized  thermo¬ 
dynamically  by  a  single  phase  (liquid)  at  temperature  T0  and  at  (applied)  pressure  P0. 

If  solidification  then  occurs  in  the  supercooled  melt,  the  system  will,  given  enough  time, 
approach  a  final  equilibrium  state  which  is  isenthalpic  to  the  initial  state.  In  other 
words,  the  assumed  conditions  of  constant  pressure  and  adiabatic  isolation  of  the  freez¬ 
ing  system  require  that  the  enthalpies  of  the  initial  and  final  equilibrium  states  be  equal. 
The  type  of  postsolidification  equilibrium  is  restricted  to  the  two  cases  predicted  from 
Gibbs’  phase  rule  for  a  unary  (one -component)  system:  (a)  invariant  equilibrium,  where 
the  final  state  is  characterized  by  two  phases  (solid  and  liquid,  at  pressure  P0),  at  their 
invariant  final  temperature  Tf  =  Tr ,  where  Tr  denotes  the  usual  equilibrium  temperature 
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between  solid  and  liquid,  and  (b)  univarlant  equilibrium,  where  the  final  state  is  charac¬ 
terized  by  one  phase  (solid,  at  pressure  P0)  at  a  temperature  Tf  <  Te,  where  here  Tf 
Is  some  function  of  the  initial  temperature  T0. 

The  relationships  between  the  thermodynamic  states  before  and  after  isenthaipic 
solidification  are  indicated  In  Fig.  1.  Figure  1  Is  a  normalized  enthalpy-temperature 
diagram  in  which  the  ordinate  is  the  usual  enthalpy  function  divided  by  the  latent  heat  of 
fusion  A0*  and  the  abscissa  is  the  absolute  temperature  divided  b>  A0/C,,,  where  C(  is 
the  specific  heat  (at  constant  pressure)  of  the  ilquid,  assumed  temperature  Independent. 
In  such  a  representation  the  (normalized)  enthalpy-temperature  relation  for  the  liquid 
phase  will  aiways  be  a  straight  line  of  unit  slope,  whereas  the  (normalized)  enthalpy  - 
temperature  relation  for  the  solid  phase  wlU  in  general  be  a  curve  with  every  point  ,is- 
piaced  vertically  from  the  ilquid  enthalpy  line  by  an  amount  A/A0.  Here,  A  is  the  latent 
heat  of  lusion,  which,  in  general,  has  a  slight  temperature  dependence.  As  indicated  in 
Fig.  1,  a  liquid  phasp  at  temperature  i?0  can  exist  in  any  of  three  types  of  contiguous, 
initial  states:  a  stable  state  (i?0  >  dj,  where  no  solidification  can  occur;  a  metastable, 
undercooled  state  (i?p  -  1  <  d0  <  de),  where  Isenthaipic  solidification  produces  invariant 
equilibrium  at  temperature  df  =  i?„;  and  a  metastable,  hypercooled  state  (d0  <  de  -  1), 
where  Isenthaipic  solidification  produces  univarlant  equilibrium  at  temperature  df  = 
i?0  +  A/A0.  The  two  metastable  states  are,  therefore,  distinguishable  on  the  basis  of 
whether  the  initial  bath  supercooling,  Ad5  tf,  -  d0,  is  less  than  unity  (undercooled)  or 
greater  than  unity  (hyper cooied).  The  effect  of  the  Initial  state  on  the  kinetics  of  freez¬ 
ing  will  now  be  considered. 


Fig.  1  -  Normalized  enthalpy- 
temperature  diagram.  Isenthai¬ 
pic  solidification  reactions  are 
indicated  by  the  initial  (■)  and 
final  (• )  states  of  the  system. 
Invariant  equilibrium  obtains  for 
the  initial  condition  de  -1  <  < 

d, ,  whereas  univariant  equilib¬ 
rium  obtains  for  the  initial  con¬ 
dition  d0  <  dr  -1.  Above  tem¬ 
perature  dp,  the  liquid  phase  is 
stable  relative  to  the  solid  phase. 


*\g  is  the  latent  heat  of  fusion  at  the  equilibrium  melting  point  i5f . 
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Kinetics, 

So  far,  oniy  the  thermodynamic  relationships  for  the  equilibrium  states  before  and 
after  isenthaipic  solidification  have  been  considered.  We  now  wiii  discuss  some  possible 
kinetic  processes  that  a  sol  ifying  system  might  undergo  as  it  departs  from  the  initial 
(metastabie)  equilibrium  state  and  approaches  a  finai  equilibrium  state. 

Planar  Solidification  at  Constant  Speed  -  Figure  2a  represents  a  system,  initially 
supercooled  an  amount  AO  =  -  60,  solidifying  under  isenthaipic  conditions  with  a 

planar  solid/iiquid  interface  advancing  at  a  constant  rate  V.  The  isenthaipic  condition 
implies  that  no  transfer  of  heat  can  occur  between  the  system  and  its  sunoundings.  if 
we  assume  that  the  thermophysicai  constants  of  the  solid  and  liquid  phases  are  independ¬ 
ent  of  temperature,  then  an  exact  solution  exists  for  the  steady-state  temperature  distri¬ 
bution  along  the  direction  of  growth,  vi2., 

0(X)  =  +  [  AO  (l  -  +  ^  ]  exp(-VX/o),  X  a  0,  (la) 

(see  Fig.  2b)  and 

0(X)  -  o0  +  [a*  (i  -  +  Xso,  (lb) 


✓  INTERFACE 


Fig.  Z  -  isenthaipic  solidification  with  a  plare- 
front  morphology:  (a)  steady-state  configuration 
of  the  system  near  the  solid/liquid  interface; 
X  =  0  denotes  interface  in  a  moving  coordinate 
system;  (b)  normalized  temperature  versus  dis¬ 
tance  in  the  vicinity  of  X  =  0 — note  that  the  melt 
supercooling  (Ad)  must  exceed  unity  as  a  neces¬ 
sary  condition  for  th's  solidification  mode;  (c) 
normalized  enthalpy  versus  distance  nearX  =  0— 
not-'  that  the  enthalpy  of  the  solid  is  equal  to  that 
of  the  liquid  far  from  the  interface.  Here,  VA„ 
is  substituted  for  the  equivalent  expression 
[Ad(l  -  Cf /Cs)  +  Cf/C,]  appearing  in  Eqs.  (la) 
and  ( lb). 
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where  X  is  the  distance  to  any  point  measured  normal  to  the  solid/liquid  interface 
(X  -  0)  on  a  coordinate  system  moving  with  the  solid  ^liquid  interface.  We  will  now  make 
the  approximation  C5  =  C(> ,  so  the  quantity  within  the  brackets  of  Eqs.  (la)  and  (lb)  be¬ 
comes  unity,  and  the  expressions  agree  with  those  given  in  Refs.  6  and  7.  The  speed  of 
advance  V  of  the  interface  is  a  constant  determined  by  the  interface  supercooling  through 
some  independent  kinetic  relationship  of  the  general  form 

V  f[*r-l«0  +  l)1-f(8).  (2) 

Here,  the  interface  supercooling  5  is  related  to  the  meit  temperature  by  the  relation 
6  =  -  (d0  +  1).  [See  Fig.  3  for  the  scale  of  (r  ormaiized)  supercooling  used  throughout 

this  report.] 


(normal  melting  point) 

(r)  (local  equilibrium 

te-.perature) 


<1?*.  (Interface 

'  temperature) 


,9  (malt 

v  °  temperature) 


Fig.  3  -  Reference  scale  of  normalized 
temperatures  and  normalized  tempera¬ 
ture  differences.  All  quantities  shown 
are  dependent  or  independent  variables 
except  i>e  which  is  a  constant. 


It  becomes  apparent  from  the  nature  of  the  heat-flow  solutions,  Eqs.  (la)  and  (lb) 
that  steady-state  advance  of  a  planar  solidification  front  under  isenthaipic  conditions  re¬ 
quires  an  initial  level  of  normalized  supercooling,  A#  =  dp  -  d0,  greater  than  unity;  i.e., 
initially,  the  meit  must  be  in  a  hypercooied  s.ate.  Moreover,  as  is  seen  in  Fig.  2c,  the 
solid  and  liquid  phases  far  ahead  of  and  behind  the  transformation  front  have  equal  en¬ 
thalpy.  Thus,  the  solidification  of  a  hypercooied  melt  (Ad  >  1)  via  steady-state  plane- 
front  transformation  appears  to  be  a  unary  analog  to  the  so-caiied  massive  transforma¬ 
tions  occurring  without  change  of  composition  in  binary  systems  (8). 

In  summary,  the  planar-fronf  solidification  just  discussed  represents  the  limiting 
case  of  isenthaipic  solidification,  wherein  the  temperature  of  every  point  on  the  solid/ 
liquid  interface  is  a  constant  amount  (unity)  above  the  meit  temperature  d0  at  infinity; 
i.e.,  =  d0  +  1.  The  interface  in  this  case  is  spatialiy  and  temporally  isothermal,  but 

the  rate  of  motion  is  determined  by  both  heat-fiow  and  kinetic  considerations. 

Nonpianar  Solidification  In  Hypercooied  Melts  -  The  instabilities  which  lead  to  the 
breakdow"  of  a  planar  solidification  front  should  be  operative  in  hypercooied  systems  as 
well  as  in  normally  um’ercooied  systems.  For  instance,  Figs.  4a  and  4b  are  illustrations 
of  nonpianar  interfaces  in  a  normally  undercooied  meit  — where  dendrites  form  — and  in 
a  hypercooied  1.  .elt  —  where  "scallops"  form.  "Scalloping"  of  the  interface  leads  to  the 
same  type  of  enhanced  i.uat  diffusion  that  dendrites  permit,  but,  with  a  scalloped  inter¬ 
face,  total  solidification  occurs  some  distance  back  from  the  most  advanced  regions  of 
the  interface.  It  will  be  shown  later  (with  light-scattering  experiments)  that  in  the  case 
of  phosphorus  the  extent  of  the  scalloped  region  diminishes  rapidly  with  increased  super¬ 
cooling.  Also,  as  the  supercooling  increases  beyond  about  Ad  >  2,  the  scalloped  inter¬ 
face  becomes,  from  a  morphological  standpoint,  almost  indistinguishable  from  a  planar 
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c)  DENDRITIC  t>)  SCALLOPED  c)  QUASI*PLANAR 

Fig.  4  -  Nonplanar  solidification  mor¬ 
phologies:  (a)  dendritic,  for  under¬ 
cooled  melts,  (b)  scalloped,  and  (c) 
quasi-planar,  for  hypercooled  melts. 

All  three  morphologies  have  complex 
thermal  diffusion  fields  associated 
with  tiieir  propagation. 


interface.  Thus,  as  the  supercooling  in  a  pure  melt  increases  beyond  unity,  a  continuous 
transition  from  dendritic,  to  scalloped,  to  quasi-planar'*'  interfaces  occurs,  c.f.,  Fig.  4. 


EXPERIMENTAL  TECHNIQUES 

White  phosphorus  (P4)  was  chosen  for  study  Decause  its  thermophysical  properties 
are  favorable  for  achieving  large  levels  of  (normalized)  supercooling.  Specifically,  a 
normalized  supercooling  of  unity  is  reached  when  liquid  P4  is  cooled  about  25  C°  below 
its  equilibrium  temperature,  Te  =  44°C;t  by  comparison,  most  metals  require  300  C°  to 
500  C°  cooling  below  their  respective  Te’s  to  reach  unit  normaliiud  supercooling.  In 
addition,  the  low  entropy  of  fusion  (9)  of  P4  makes  it  likely  that  some  aspects  of  its 
solidification  behavior  are  similar  to  those  of  the  pure  metals  (10). 


Establishment  of  the  Transition  Temperature 

A  detailed  description  of  the  apparatus  a  id  experimental  techniques  used  in  deter¬ 
mining  the  transition  temperature  between  invariant  and  univuriant  isenthalpic  solidifica¬ 
tion  in  P4  has  been  given  elsewhere  (11).  It  will  serve  the  purpose  of  this  report,  how¬ 
ever,  to  repeat  some  of  those  results  whicu  are  pertinent  here.  Figure  5  is  a  plot  of 
measured  postrecalescence  temperature  Tf  versus  Initial  molt  temperature  T0  for  a 
specimen  of  vacuum  distilled  P4  undergoing  isenthalpic  solidification.  The  data  shown 
in  Fig.  5  exhibit  a  sharp  slope  discontinuity  at  T0  =  16.4° C.  In  the  temperature  range 
18.4°C  <  T0  <  44°C  the  solidification  was  of  the  normal  invariant  type;  i.e.,  the  post- 
recalescence  temperature  was  independent  of  the  initial  temperature  and  was  always 
equal  to  the  (invariant)  equilibrium  temperature  Te.  However,  for  T0  <  18.4°C  the 
postrecalescence  temperature  failed  to  reach  Te  —  falling  ihort  about  1.06  C°  for  each 
additional  degree  of  initial  supercooling.  This  latter  behavior  is,  of  course,  characteris¬ 
tic  of  univariant  solidification,  where  Tf  is  a  function  of  T0.  This  experiment  served  to 


''The  term  quasi-planar  denotes  macroscopic  planarity  of  the  interface,  but  on  the  mi¬ 
croscopic  scale  of  the  heat-ilow  the  interface  does  not  necessarily  behave  as  a  plane 
front. 

tThe  symbol  C°  indicates  a  temperature  interval,  and  the  symbol  °C  indicates  a  specific 
temperature. 
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Fig.  5  -  Final  temperature  T  ,  established  after 
isenthalpic  solidification  of  'supercooled  liquid 
phosphorus,  versus  initial  temperature  T„  of  the 
molten  phosphorus.  The  equilibrium  temperature 
Tt  is  44°C,  and,  as  denoted  by  the  sudden  „  angc 
in  slope,  the  critical  temperature  for  univariant 
solidification  T,  -  X0/Cf  is  '.  S .4  *  0.1°C. 


establish  unequivocally  the  minimum  level  of  supercooling-  required  for  univariant  isen¬ 
thalpic  solidification  as  ATmin  =  25.6  C°. 


Light  Severing  During  ,  olidificat<on 

An  apparatus  was  constructed  to  permit  an  evaluation  of  structural  changes  during 
solidification  in  both  the  undercooled  and  hypercooled  temperature  ranges,  viz.,  T0  above 
and  below  18.4°C,  respectively.  This  apparatus  consisted  of  a  light-tight  structure  con¬ 
taining  a  stabilized  collimated  light  source,  a  highly  sensitive  rapid-respor.se  CdS  photo¬ 
detector  and  readout  circuitry,  a  specimen  chamber,  and  a  "trigger"  mechanism  for  in¬ 
ducing  solidification  In  the  P4  specimen  at  any  desired  level  of  supercooling.  An 
experiment  consisted  of  measuring  the  light  intensity  transmitted  by  a  specimen  of 
supercooled  liquid  P4  and  then  monitoring  the  changes  in  transmitted  intensity  as  freez¬ 
ing  took  place.  The  change  in  transmitted  light  intensity  was  caused  by  scattering  of  the 
mala  beam  on  solid/liquid  and  solid/solid  interlaces  that  formed  within  the  specimen. 
Thus,  the  decrease  in  transmitted  light  intensity  was  related  to  the  average  density  of 
internal  interlaces  along  the  light  path.  Figure  6  is  a  composite  diagram  of  scattered 
light  intensity  versus  time,  for  various  initial  temperatures.  Just  prior  to  "triggering" 
solidification  In  each  experiment  (t  =  0  min),  the  transmitted  light  intensity  was  adjusted 
to  an  arbitrary  scale  of  100,  or,  equivalently,  the  scattered  light  intensity  was  set  to 
zero  on  an  arbitrary  scale  of  100.  Thus,  all  the  curves  of  Fig.  6  have  common  ordinate 
and  abscissa  axes.  Because  the  incident  beam  intensity  war  held  constant  during  solidi¬ 
fication,  it  is  most  convenient  to  think  of  the  ordinate  scale  in  Fig.  6  as  proportional  to 
the  density  of  scat* Ing  sites  present  during  freezing. 
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Fig.  6  -  Percent  scattered  light  in¬ 
tensity  versus  time  after  initiation 
of  solidification  in  supercooled  P4 . 
Curves  A,  B,  C,  and  D  are  for  un¬ 
dercooled  melts  with  progressively 
increasing  supercooling.  Curves  E, 
F,  and  G  are  for  hypercooled  melts 
with  progressively  increasing  super¬ 
cooling.  Note  the  rapid  change  in  be¬ 
havior  as  the  initial  temperature 
drops  beloe  the  critical  temperature 
for  univariant  isenthalpic  solidifica¬ 
tion  ( 18.4°C). 
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We  see,  starting  with  the  ieast  underoooled  (warmest)  melts  (curves  /’  and  B),  that 
freezing  caused  a  rapid  increase  in  scattering  sites,  followed,  thereafter,  by  o-  iy  slight 
and  somewhat  irreguiar  changes.  This  is  typical  of  normal  dendritic  freezing,  where  the 
finely  dispersed  solid/liquid  interfaces  formed  du  .'ing  passage  of  the  main  solidification 
front  persist  for  relatively  iong  times.  If  the  system  were  truly  adiabatic,  the  solid/ 
liquid  mixture  would,  in  time,  coarsen  slightly,  and  the  scattered  intensity  would  diminish 
very  siowiv.  Some  heat  losses  and  further  solidification,  however,  were  unavoidable  over 
the  period  of  the  experiment  (usuaiiy  3  to  5  min),  and  the  decrease  in  scattering  sites 
(solid/iiquid  interface  area)  accompanying  this  additional  (nor.isenthalpic)  solidification 
accounts  for  the  behavior  of  these  curves.  For  specimens  frozen  at  19.5°C  and,  even 
more  pronounced,  for  those  frozen  at  18.6°C,  the  scattering-site  density  shows  a  con¬ 
spicuous  maximum  followed  by  a  rapid  decrease  to  a  nearly  constant  level.  This  behav¬ 
ior  is  characteristic  of  systems  solidifying  from  liquids  with  normalized  supercooling 
near  unity,  where  the  enthalpy  of  the  liquid  phase  at  the  initial  temperature  is  approxi¬ 
mately  equal  to  the  finai  enthaipy  of  the  solid  phase  at  the  equilibrium  temperature.  Just 
after  passage  of  the  main  solidification  front,  a  transitory  solid  and  liquid  structure  ex¬ 
ists  containing  numerous  pockets  of  interdendritic  liquid;  subsequently,  *hese  liquid  re¬ 
gions  solidify  more  or  Jess  rapidly  by  transferring  their  heat  of  fusion  locally  to  the 
solid  pJiase.  This  secondary  solidification  process  leads  to  a  dramatic  decrease  in  the 
number  oi  scs  ttering  sites,  and  either  goes  to  completion  (if  the  initial  supercooling  was 
siightly  greater  than  unity)  or  stops  when  the  system  finally  reaches  the  equilibrium 
temperature  (if  the  initial  supercooling  was  slightly  less  than  unity). 

As  the  initial  normalized  supercooling  of  the  specimen  rises  above  unity,  that  is,  as 
hypercooling  begins,  the  maximum  in  the  scattered  light  intensity  diminishes  in  ampli¬ 
tude  and  width,  untii  in  the  case  of  the  curves  for  temperatures  below  16.9°C,  the  maxima 
are  no  ionger  resolved  on  the  scale  used  in  Fig.  6.  The  solidification  structures  respon¬ 
sible  for  the  lower  curves  (F,  G)  are  qualitatively  similar  to  those  responsible  for  curves 
D  and  E,  except  that  no  interdendritic  liquid  persists  in  the  structures.  Instead,  it  ap¬ 
pears  that  the  dendritic  nature  of  the  solidification  must  be  confined  to  a  fairly  narrow 
zone  just  behind  and  moving  with  the  main  solidification  front.  This  last  case  of  freezing 
in  a  pure  hypercooled  liquid  corresponds  to  what  was  previously  termed  a  scalloped  fr  • 
terface  (Fig.  4). 
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Speed  of  Solidification 

Figure  7  is  a  diagram  of  the  apparatus  employed  in  measuring  the  linear  rates  of 
solidification  in  supercooled  liquid  P4.  Two  tungsten-Kovai  ■>  thermocouples  (W1/K1 
and  W2/K2)  were  used  to  determine  the  specimen's  temperature  prior  to  freezing  and  to 

verify  that  no  appreciable  temperature  gradients 
existed  along  the  10-cm  growth  path.  To  meas- 
ure  the  time  for  the  solidification  front  to  tra- 
^  V  verse  the  growth  path,  the  two  thermocouples 

Ci0j  A f  were  connected  in  series  (opposing  polarity)  as  a 

—  solution  I  I  differential  thermocouple,  the  output  signal  from 

ft  52  \  /_\_2  which,  after  suitable  preamplification,  was  dis- 

HZ. rxF  played  on  an  oscilloscope. 


-  PYREX 
CHAMBER 


Fig.  7  -  Apparatus  for  deter¬ 
mining  linear  solidification  The  average  speed  of  solidification  was  cal- 

rates  in  supercooled  phospho-  a6  the  length  of  the  growth 

specimen  can  be  cleansed  by  f»th  to  the  transit  time  measured  from  the  wave- 

periodically  flushing  Cr03  form  produced  by  the  thermocouples.  Figure  9 

solution  through  the  apparatus.  shows  a  logarithmic  plot  of  the  speed  of  solidifi¬ 

cation  versus  the  degree  of  normalized  super¬ 
cooling  prior  to  freezing.  The  most  significant 
feature  of  these  results  is  that  no  discontinuity  in 
slope,  or  other  peculiarity,  occurred  in  the  velocity-supercooling  relationship  as  the  sys¬ 
tem  became  hypercooled,  i.e.,  as  -  #0  exceeded  unity.  It  willdBPbe  shown  that  inter¬ 
face  temperatures  can  be  estimated  when  solidification  rates  are  known  in  the  hypsrcooled 
temperature  range,  but  considerable  analysis  is  required  to  achieve  this  end. 


A  rate  measurement  was  begun  in  the  uni¬ 
formly  supercooled  specimen  by  flushing  a  stream 
of  cold  helium  gas  through  the  Kovar  tube  that 
passed  through  the  base  of  the  15-mm-diameter 
Pyrex  column.  The  resulting  rapid  cooling 
chilled  the  P4  specimen  locally  to  the  nucleation 
temperature;  moreover,  the  chilling  was  applied 
so  suddenly  that  no  significant  distortion  was  in¬ 
troduced  into  the  uniform  temperature  field  be¬ 
tween  the  thermocouple  junctions.  The  P4  crys¬ 
tals  that  nucleated  around  the  chilled  Kovar  tube 
spread  rapidly  through  the  supercooled  specimen, 
and  the  abrupt  change  in  temperature  accompany¬ 
ing  the  solidification  front  was  sensed,  in  turn,  by 
junctions  W2/K2  and  W1/K1.  Figures  8a  and  8b 
are  representative  waveforms  of  the  signal  gen¬ 
erated  by  the  passage  of  the  solidification  front 
across  the  junctions  of  the  differential  thermo¬ 
couple.  The  fast  response  achieved  with  this  sys¬ 
tem  (detectable  signals  within  1  millisecond)  is 
attributed  primarily  to  the  small  gauge  of  the 
tungsten  and  Kovar-A  wires  (0.005-in.  diameter) 
and  to  the  fact  that  since  P4  is  an  electrical  insu¬ 
lator,  no  additional  insulation  was  needed  on  the 
junctions. 
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Fig.  8  -  Representative  waveforms  of  the  output  signal 
from  the  solidification-rate  apparatus.  The  transit 
time  for  interface  advance  between  the  thermocouple 
junctions  it  taken  as  the  interval  between  the  "trig¬ 
gering"  point  of  the  waveform  and  the  abrupt  {negative) 
increase  in  slope.  Transit  times  in  phosphorus  varied 
between  25  msec  and  1  sec  for  the  range  of  supercool¬ 
ing  studied 


o)  Mslt  Umpsraturs :  34.4  *C 
V  •  43.4  cm  /  tee 


10  IBM#  /4I*  1  » 

b)  Mslt  tsmpsralurt1  14.4  *C 
V*  204  cm/M* 


M?T  TEMPERATURE  [*C] 


Fig.  9  -  Velocity  of  the  solidification  front  versus 
melt  so  ercooling  (lower  abscissa)  and  melt  tem¬ 
perature  (upper  abscissa).  The  curvature  of  the 
plot  un  a  log-log  scale  indicates  that  no  simple 
power-law  relationship  exists  between  V  and  Ad 
over  the  range  of  Ad  investigated  here. 
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DATA  ANALYSIS 
Planar  Growth  Hypothesis 

Horvay  and  Cahn  (12)  have  shown  In  a  rigorous  fashion  that  when  a  crystal  grows 
into  a  uniformly  supercooled  melt,  the  maximum  normalized  temperature  difference 
that  can  exist  between  the  solid/liquid  interface  (assumed  spatially  isothermal)  and  its 
melt  (far  from  the  interface)  is  unity.  Moreover,  irrespective  of  interface  morphology, 
if  this  normalized  temperature  difference,  dt  -  d0  (see  Fig.  3)  approaches  unity,  then 
interface  motion  and  shape  become  characterized  by  a  P£clet  number  that  approaches 
infinity.  The  Peclet  number  Pe  is  a  convenient  dimensionless  heat-flow  parameter  de¬ 
fined  by  the  relationship 


where  V  and  r  are  the  normal  velocity  and  mean  radius  of  curvature  of  the  interface, 
respectively,  and  a  is  the  thermal  diffusivity  of  the  melt.  V,  of  course,  must  always  be 
finite,  so  a  Peclet  number  approaching  infinity  implies  that  r  is  approaching  infinity,  or, 
equivalently,  that  the  mean  gaussian  curvature  of  the  Interface  is  approaching  zero.  This 
general  conclusion  reached  by  Horvay  and  Cahn  is  formally  equivalent  to  the  specific  re¬ 
sult  shown  by  Eqs.  (la)  and  (lb)  In  our  discussion  of  planar  isenthalpic  solidification: 
viz.,  that  a  planar  solidification  front  propagates  through  the  melt  with  a  unitary  normal¬ 
ized  temperature  rise  above  the  melt  temperature,  and  that  the  planar  solidification  front 
is  the  most  sluggish  type  of  front  Insofar  as  a  maximum  amount  of  the  available  super¬ 
cooling  is  used  up  in  driving  the  heat  transfer  process.  Since  in  our  experiments  we 
were  able  to  achieve  normalized  supercooling  levels  as  large  as  1.7,  the  possibility  must 
be  considered  that  planar  solidilication  could  take  place  in  ?4 ,  with  the  true  interface 
temperature  as  much  a-,  18.4  C°  (0.7  normalized)  below  the  equilibrium  temperature. 

In  planar  solidification,  where  no  Gibbs-Thompson  effect  occurs,  the  normalized 
supercooling  6  responsible  for  kineti:  processes  (other  than  heat  flow)  is  given  by 

6  =  de  -  =  Ad  -  1.  (3) 

A  logarithmic  plot  o*  solidification  speed  V  against  Ad  -  1  would  either  reveal  the  true 
dependence  of  V  on  6  (if  the  morphology  was  Indeed  microscopically  planar)  or  prov-' 
unequivocally,  that  the  morphology  was  nonplanar  even  at  extreme  levels  of  supercooling. 
Theoretical  works  by  Wilson  (13),  Frenkel  (14),  Hillig  and  Turnbull  (15),  and  more  re¬ 
cently,  by  Cahn  (16)  and  by  Cahn,  Hillig,  and  Sears  (17)  have  shown  the  existence  of  either 
a  quadratic  or  linear  relationship  between  V  and  5,  depending  on  whether  interface  mo¬ 
tion  occurred  by  lateral  spreading  of  growth  layers  from  dislocation  sources,  or  by  uni¬ 
form  advancement  via  random  atomic  attachment,  respectively.  Figure  10  indicates  that 
neither  a  linear,  quadratic,  nor  other  simple  power-law  relation  exists  between  V  and 
Ad  -  1.  Instead,  we  conclude  from  these  data  that  6  *  Ad  -  1,  and,  therefore,  that  P4 
solidifies  with  a  nonplanar  (microscopically  scalloped)  interface  even  when  the  total 
supercooling  is  as  large  as  1.7.  The  presence  of  a  nonplanar  interface  during  solidifica¬ 
tion  demands  that  both  the  Gibbs-Thompson  effect  and  the  effect  of  multidimensional  heat 
flow  be  considered  when  attempting  to  deduce  tl  interface  kinetic  relationship  [V  =  f(6)l 
from  the  experimental  rate  data  [V  =  g(Ad)]. 


Analysis  for  Nonplanar  Growth 

In  the  preceding  section  re  proved  indirectly  that  P4  freezes  with  a  nonplanar  inter¬ 
face  morphology  from  both  undercooled  a.;d  hypercooled  melts.  To  proceed  further  wich 
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Fig.  10  -  Growth  velocity  (corrected  for  viscosity) 
a-  a  function  of  interface  temperature  Ai>  -  1,  as¬ 
suming,  a  priori,  a  planar  interface  morphology 
during  growth  from  hypercooled  melts.  The  ab¬ 
sence  of  a  linear  or  quadratic  relationship  between 
V,.  and  A 0  -  1  indicates  that  no  planar  growth  actu¬ 
ally  occurs. 


an  interface  temperature  determination,  some  estimate  must  be  made  of  the  shape  of  the 
solid/liquid  interface  in  these  regions*  where  the  interface  advances  normal  to  itself 
most  rapidly. 

Phosphorus  (a  phase)  is  composed  ot  rather  symmetrical  P4  molecules  arranged  on 
a  bcc  lattice  (18)  with  a  large  unit  cell  that  contains  56  P4  molecules  (19).  The  rotations 
of  the  P.  molecule  remain  unquenched  in  the  solid  phase  down  to  low  temperatures 
(196°K)  (20,21,  which  accounts  for  phosphorus'  low,  typically  metallic,  entropy  of  fusion 
(1.98  cal/°K-mole  P4).  In  view  of  these  properties,  we  believe  that  in  P4  the  protuber¬ 
ances  (dendrites  or  scallops)  that  constitute  the  nonpl  -»nar  solid/liquid  interface  are 
similar  to  those  found  in  metals  (2).  Moreover,  we  shall  assume,  a  priori,  that  a  parab¬ 
oloidal  surface  adequately  approximates  the  true  shape  of  the  pro.  iberances  over  a  lim¬ 
ited  region  near  their  tips. 

It  has  been  shown  by  Bolling  and  Tiller  (22)  that  an  isothermal  paraboloid  cannot  be 
a  precise  description  of  the  state  of  an  isolated  growing  dendrite;  the  sides  of  the  den¬ 
drite  have  leso  curvature  than  the  tip  and  grow  more  slowly,  and  must,  therefore,  be 
warmer  than  die  tip.  However,  when  we  consider  g  owth  at  large  normalized  supercool¬ 
ings,  adjacent  protuberances  will  tend  to  warm  the  sides  of  one  another  and  thus  make 
more  probable  the  propagation  of  approximately  paraboloidal  shapes  with  the  required 
nonisothermal  surface.  In  addition,  Table  i  shows  that  for  P4  the  ratio  of  the  thermal 
diffusion  length  a/V  to  the  critical  radius  for  nucleation  r*  is  only  about  14,  over  most 


*The  tip  area  of  a  dendrite  or  scallop  is  an  example  of  these  regions;  in  particular,  note 
that  measurement  of  solidification  speed  usually  is  confined  to  measurement  only  of  tip 
speed. 
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of  the  range  of  experimentally  attainable  supercooling.  Since  r*  is  the  minimum  radius 
possible  for  a  dendrite  tip,  the  small  a/Vr*  values  indicate  that  the  tip  region  of  any 
particular  interface  protuberance  has  a  thermal  environment  that  is  not  strongly  affected 
by  the  thermal  diffusion  fields  of  similar  neighboring  protuberances,  nor  is  it  affected 
significantly  by  any  shape  instabilities  developing  behind  the  tip  region  (4).  F"U  these 
reasons,  a  tractable,  yet  not  unreasonable,  heat-flow  analysis  can  be  applied  here, 
wherein  the  tip  regions  of  the  interface  protuberances  arc  considered  as  independently 
advancing  paraboloids,  the  tip  temperatures  of  which  are  not  affected  by  the  heat  released 
by  neighboring  protuberances. 


Table  1 

Calculated  Quantities  for  P4 


Super¬ 

cooling 

Ad 

Critical 
Radius  r* 

(cm  x  107) 

Thermal 
Diffusion 
Length  a/V 
(cm  x  10s) 

a 

Vr* 

0.2 

27.5 

12.3 

45 

0.4 

13.8 

3.2 

23 

0.6 

9.2 

1.5 

16 

0.8 

6.9 

0.91 

13 

1.0 

5.5 

0.78 

14 

1.2 

4.6 

0.63 

14 

1.4 

3.9 

0.53 

14 

1.6 

3.4 

0.47 

14 

.  1.7 

3.2 

0.46 

14 

Values  of  solidification  speed  were  read  from  the  smooth  curve  drawn  through  the 
V  vs  Ad  data  in  Fig.  9,  at  normalized  supercooling  values  that  conveniently  covered  the 
range  of  observation,  viz.,  0.2,  0.4,  0.6,  0.8,  1.0,  1.2,  1.4,  1.6  and  1.7.  Each  coordinate 
pair  (V,  Ad)  formed  the  basis  for  a  diagram  such  as  shown  in  Fig.  11.  Figures  11a  and 
lib  relate  normalized  temperature  differences  to  the  tip  radius  of  the  paraboloidal  pro¬ 
tuberances  which  we  assume  describe  the  tip  regions  of  the  solid/liquid  interface.  The 
tip  radii  are,  for  convenience,  expressed  in  units  of  the  critical  radius  for  nucleation.* 
The  curves  labeled  GT,  representing  the  difference  between  tip  equilibrium  temperature 
de(r)  and  melt  temperature  d0  at  infinity  were  calculated  as  a  function  of  r/r*  v/ith  the 
Gibbs -Thompson  equation, 


*e(r)  -  *o  =  -  r*/r); 

see  Fig.  3  for  the  definitions  of  each  normalized  temperature  or  temperature  difference. 
The  curves  labeled  HC  represent  the  differences  between  normalized  interface  tempera¬ 
ture  d j  and  .nelt  temperature  d0  at  infinity  calculated  from  the  Horvay-Cahn  heat-flow 
equation  (12)  for  circular  paraboloids, 

dj  -  d0  =  -P£[exp(Pe)]  Ei(-Pe).  (4) 

Although  the  temperature  difference  given  by  Eq.  (4)  is  most  conveniently  expressed  as  a 
function  of  the  P^clet  number  P4  it  should  be  noted  that  the  variables  Px  and  r/r*  are 
simply  related  through  the  identity 


*Critical  vadii  ere  calculated  here  on  the  basis  of  a  solid/liouid  surface  energy  of 
8  erg/cm  J. 
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Fig.  11  -  Temperature  differences  aris¬ 
ing  from  curvature  and  heat-flow  effects 
at  the  tips  of  advancing  paraboloids,  as 
functions  of  the  tip  radius,  (a)  Under¬ 
cooled  melt  (Afl  =  0.6),  showing  finite 
range  of  acceptable  values  of  r/r*  which 
give  positive  values  of  5.  (b)  Hyper- 

cooled  melt  (Ad  =  1.6),  showing  a  semi¬ 
infinite  range  of  acceptable  values  of 
r/r*.  Ad  represents  the  position  of  the 
normal  melting  point  relative  to  the  tem¬ 
perature  of  the  supercooled  bath;  the 
curves  designated  GT  (Gibbs-Thompson) 
represent  the  local  equilibrium  temper¬ 
ature  de(r)  as  affected  by  the  tip  curva¬ 
ture;  the  curves  designated  HC  (Horvay- 
Cahn)  represent  the  rise  of  interface 
temperature  uk  above  melt  temperature, 
due  to  latent  heat  emission.  The  HC 
curves  are  based  on  the  measured  s;rowth 
velocity  observed  in  P4  for  each  value 
of  Ad.  5  is  the  interfacial  supercooling 
driving  solidification  kinetics. 


— 

2a  \r */ 


where  the  specific  constant  of  proportionality  between  Pe  and  r/r*  is  again  derived 
from  the  experimental  values  of  V  and  Ad  upon  which  the  particular  diagram  is  based. 
Hence,  each  HC  curve  expresses  d ;  -  d0  as  a  function  of  r/r* . 

Since  the  ordinate  difference  between  curves  GT  and  HC  represents  exactly  the  nor¬ 
malized  temperature  difference  6  responsible  for  "driving,"  at  the  tip,  kinetic  processes 
other  than  heat  flow,  it  is  clear  then  that  Figs,  11a  and  lib  succinctly  reveal  the  total 
range  of  solidification  behavior  consistent  with  our  model  and  with  experiment.  In  other 
words,  for  each  coordinate  pair  (V,  Ad)  a  specific  range  of  tip  radii  is  obtained  over 
which  solidification  is  possible,  that  is,  over  which  nonzero,  positive  values  exist  for  6. 

A  striking  difierence  between  the  behavior  of  a  normally  undercooled  melt  (Fig.  11a)  and 
a  hypercooled  melt  (Fig.  lib)  is  that  the  solidification  must  occur  over  a  limited  .  ange  of 
r/r*  for  the  former,  whereas  a  semi-infinite  range  of  r/r*  is  admissible  for  the  latter. 


Determination  of  Interface  Temperature 

Cahn,  Hillig,  and  Sears  have  shown  that  the  kinetics  of  atomic  attachment  at  a  solid/ 
liquid  interface  can  be  displayed  graphically  by  plotting  V[»j(d,  )/tj0]/6  versus  6,  where 
the  factor  t]('Si)/T)Q  is  the  ratio  of  the  melt's  viscosity  at  the  interface  temperature  d;  to 
the  melt's  viscosity  at  the  normal  melting  temperature.  Correction  of  the  velocity  by 
the  ratio  ^(d()/?;0  compensates  for  the  changes  in  the  intrinsic  solidification  rate  caused 
by  changes  in  the  self-diffusion  rate  with  temperature  in  the  region  of  the  solid/liquid 


> 
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interface.  Of  course,  to  apply  their  analysis  we  first  must  determine  the  particular 
value  of  or  r/r*  (since  6  i?  already  a  known  function  of  r/r*)  which  characterizes 
the  P4  system  at  each  (v,  Ail)  coordinate. 

Figure  12  is  a  composite  plot  of  V[tXi>i)/??0]/6,  henceforth  denoted  Vc/6,  versus 
r/r*  for  various  fixed  values  of  the  experimental  coordinates  Ad  and  V.  Tl»e  self¬ 
diffusion  corrections  used  in  Fig.  12  are  based  on  published  values  for  the  viscosity  of 
P4  (23,24).  These  corrections  were  calculated  on  the  basis  of  the  interface  temperatures 
d  j  computed  by  the  Horvay-Cahn  analysis  for  each  value  of  r/r*,  Ad,  and  V.  The  curves 
in  Fig.  12  represent  the  possible  combinations  of  Vc  /6  and  r/r*  which  are  consistent 
with  the  experimental  measurements  of  V  and  Ad,  and  with  the  paraboioidal  model  of 
the  solid/liquid  interface.  We  see  that  the  minimum  values  of  Vc/6  all  fall  within  a  few 
percent  of  each  other;  that  is,  (Vc/6)m;n  is  approximately  constant  over  a  wide  range  of 
Ad.  In  addition,  the  curves  drawn  for  the  hypercooled  meits  indicate  that  .he  maximum* 
values  of  Vc/6  are  decreasing  with  increased  supercooling  and,  moreover,  are  approach¬ 
ing  (vc/6)min- 

Strictly  speaking,  the  foregoing  analysis  provided  only  upper  and  lower  limits  for 
Vc/6  over  the  entire  supercooled  range.  Fortunately,  in  P4  we  have  found  that  the  upper 
limit  converged  rapidly  toward  a  neariy  constant  iower  limit  as  Ad  increased.  We  con¬ 
clude  from  this  behavior  that  the  values  of  (Vc/6)min  are  the  characteristic  values  we 
seek  to  describe  the  solidification  behavior  of  this  system,  and  that  the  corresponding 
values  of  r/r*  represent  the  true  tip  radii  of  the  dendrites  or  scallops.  A  satisfying  as¬ 
pect  of  this  result  is  that  without  recourse  to  any  a  priori  maximization  principle  we 
have  found  strong  indications  that  P4  solidifies  with  the  maximum  amount  of  interfaciai 
supercooling  consistent  with  the  stipulated  model  for  interface  morphology.  Solidification 
at  maximum  kinetic  supercooling  with  a  given  Ail  is  formally  equivalent  to  solidification 
at  maximum  veiicity  with  a  given  Vc  /6.  However,  since  no  legitimate  theoretical  basis 
exists  for  the  maximum  velocity  principle  (25),  these  results  for  P4  are  given  only  as 
empirical  findings  and  not  as  a  justification  for  the  general  applicability  of  any  maximi¬ 
zation  principle.  Obviously,  this  point  requires  additional  theoretical  review  as  well  as 
further  experimentation  on  nther  systems. 


Interface  Kinetics  for  P4 

If  we  accept,  for  the  tine,  (Vc/6)min  as  the  characteristic  parameter  for  solidifica¬ 
tion  in  P4,  then  the  kinetic  behavior  can  be  represented  conveniently  by  a  plot  of  (Vc/6)min 
against  6,  as  shown  in  Fig.  13.  The  upper  curve  in  Fig.  13  showing  (Vc/6)max  versus  6 
was  added  to  emphasize  the  rapid  convex  of  this  upper  bound  toward  (Vc/6)mjn  as 
6  increases.  It  is  apparent  from  Fig.  13  t*.*t  a  iinear  dependence  exists  between  Vc  and 
6,  since  (Vc/6)min  is  a  constant  (above  6  =  0.1).  The  vaiue  of  this  constant  is  450  ±  10 
cm/sec,  corresponding  to  a  rate  constant  of  17.7  ±  0.4  cm/sec  C°.  The  iinear  relation¬ 
ship  between  Vc  and  6  implies  that  uniform  attachment  kinetics  are  controiiing  at  the 
soiid/liquid  interface  in  P4  when  the  meit  supercooling  Ad  is  greater  than  about  0.35 
(or  on  a  dimensional  scaie,  greater  than  about  9  C°).  The  data  in  Fig.  13  also  show  some 
indication  of  the  onset  of  the  so-caiied  transitional  growth  regime  beiow  about  6  =  0.1. 
Below  6  =  0.1,  the  curve  bends  downward  rather  abruptly  as  the  growth  mechanism  pre¬ 
sumably  changes  from  uniform  attachment  kinetics  (linear  in  5)  to  iayer-passage  limited 
kinetics  (quadratic  in  6). 


*We  specifically  exclude  from  the  present  analysis  the  steeply  rising  values  ofVc/S  at 
values  of  r/r*  near  unity.  Steady-state  solidification  at  small  r/r*  is  considered  phys¬ 
ically  untenable,  insofar  as  the  kinetic  supercooling,  5,  would  be  very  unstable  relative 
to  any  slight  variations  ir.  interface  curvature,  which  occur  during  solidification.  Steady- 
state  solidification  requires  a  reasonably  weak  dependence  of  Vc  /i  on  r/r*;  for  this 
reason,  we  do  not  believe  that  steady-state  growth  can  occur  with  r/r*  near  unity. 
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Fig.  12  -  Curves  showing  Vc/8  as  a 
function  of  r/r*,  for  several  values 
of  the  melt  supercooling  Ad.  8  is  cal¬ 
culated  as  in  Fig.  1 1  by  using  the 
measured  value  of  V  corresponding  to 
each  Ad.  Curves  for  Ad  >  0.4  all  show 
minimum  values  of  Vc/8  at  about  450 
cm/sec.  In  addition,  curves  for  Ad  >  1 
yield  maximum  values  for  Vc/S  at 
large  r/r*. 


t 


Fig,  13  -  Minimum  and  maximum  values  of 
Vc /8  ,  from  curves  in  Fig.  12,  shown  as  func¬ 
tions  of  8.  Note  that,  as  8  increases,  the  val¬ 
ues  of  (Vc  /8  )„  „  appear  to  approach  the  almost 
constant  level  of  (Vc/8)mtn.  The  convergence 
of  (Vc/&)wt„  and  (Vc /S  )ltI  limits  the  uncer¬ 
tainty  in  the  true  value  of  (Vc/8)  to  less  than 
30  percent.  However,  it  is  believed  that 
(V./S)  is,  in  fact,  the  true  value. 

c  min 
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Our  preliminary  observations  of  solidification  in  P4  at  melt  supercooling  near  1  C 
confirm  the  presence  of  faceted  growth  morphologies*  at  very  iow  ieveis  of  supercooling. 
It  seems  clear,  therefore,  that  somewhere  in  the  range  of  meit  supercooling  between  1  C° 
and  9  C°  the  atomic  mechanism  of  crystal  growth  in  P4  changes  from  iayer  spreadi  r  to 
uniform  attachment.  The  present  series  of  experiments  are  qualitatively  explained  by 
the  molecular  crystal  growth  theory  of  Cahn,  Kiiiig,  and  Sears,  and  are  not  understood  in 
terms  of  Jackson's  theory  (26)  of  smocth  and  rough  interfaces  —  a  theory  which  does  not 
predict  the  observed  change  in  growth  mechanism  with  increased  kinetic  driving  force. 

In  their  paper,  Cahn,  Hillig,  and  Sears  accepted  the  approximately  quadratic  relation¬ 
ship  between  Vc  and  Ad,  which  had  been  observed  in  P4  by  Poweii,  Giiman,  and  Hiide- 
brand  (27)  as  evidence  for  layer-passage  limited  (nonlinear)  kinetics.  Although  the  pres¬ 
ent  experiments  are  in  agreement  with  Poweii,  Gilman,  and  Hildebrand's  observations, 
we  see  no  merit  in  Cahn,  Hiliig,  and  Seals'  justify ...ation  that  this  quadratic  relationship 
indicates  nonlinear  kinetics,  or,  equivalently,  tha.  the  melt  supercooling  Ad  is  nearly 
equal  to  the  true  interface  supe.  cooiing  6.  Our  analysis  has  shown,  instead,  that  6  and 
Ad  are  related  by  rather  involved  formulas,  and,  as  indicated  in  Fig.  14,  6  and  Ad  dif¬ 
fer  appreciably  in  magnitude  in  P4 .  The  'complicated  dependence  of  6  and  Ad  indicated 
in  Fig.  14  S'  iffices  to  undermine  any  straightforward  interpretation  of  the  relationship 
between  Vc  and  Ad  except,  perhaps,  at  very  low  values  of  Ad  where,  for  faceted  growth, 
the  interface  temperature  might  be  approximated  reasonably  by  the  ambient  temperature 
of  the  melt.  Fortunately,  the  relationship  determined  between  Vc  and  6  is  interpreted 


MELT  SUPERCOOLING  (4$) 

Fig.  14  -  Relationship  of  interface  super¬ 
cooling  i  [derived  from  (Vc/i>imin]  to  total 
melt  supercooling  The  plot  indicates 

that  &  is  considerably  less  than  /\d  and  that 
r.o  simple' relation  exists  between  them. 


*Faceted  growth  morphologies  are  strong  evidence  for  layer-passage  limited  kinetics. 
See  Ref.  17  for  additional  information  on  this  point. 
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easily,  except  where  transitional  growth  occurs  at  low  values  of  6.  More  work  with  P4 
is  needed  in  the  lower  ranges  of  supercooling  over  wnich  the  atomic  mechanism  for 
crystal  growth  changes  character. 

Figure  15  shows  the  fractional  distribution  of  total  (melt)  supercooling  into  its  three 
components,  6,  dt  -  d0,  and  d„  -  dc(r).  A  surprising  feature  disclosed  here  is  that  the 
fractional  supercooling  for  interface  attachment  processes  6/ Ad  rises  rapidly  with  A# 
and  accounts  for  more  than  half  the  total  supercooling  above  Ail  =  0.9.  Even  at  the  lower 
levels  of  supercooling  for  dendritic  growth,  6  accounts  for  more  than  10  percent  of  the 
total. 


Fig.  15  -  Fractional  distribution  of 
melt  suDercooling  Ad  ascribed  to  ki¬ 
netics  (i>/Ad],  heat  flow  [(iJ|  -  d0)/Ad], 
and  curvature  [(d„  -  de(r))/Ad],  as 
functions  of  Ad.  Note  the  greatly  in¬ 
creased  importance  of  the  kinetic 
contribution  at  large  supercooling. 
The  curves  should  not  be  extrapolated 
to  supercoolings  much  below  Ad  =  0.2, 
because  faceted  growth  becomes  prev¬ 
alent  at  low  supercooling. 


The  fractional  supercooling  required  for  heat  flow  (di  -  d0)/Ad  accounts  for  about 
half  the  total  supercooling  at  low  values  of  Ad,  but  only  for  about  25  percent  at  he  larger 
values  of  Ad.  Similarly,  the  fractional  supercooling  "lost"  because  of  curvature  effects 
[de  -  d„(r)]/Ad  decreases  from  over  30  percent  at  small  supercooling  to  a  nearly  con¬ 
stant  10  percent  at  large  supercooling. 

From  the  data  in  Fig.  15  we  conclude  that  heat  flow,  attachment  processes,  and  cur¬ 
vature  effects  individually  represent  important  phenomena  which  collectively  determine 
the  overall  solidification  process  In  P4;  no  one  effect  is  extraordinarily  dominant  or 
rate  limiting. 
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The  thermodynamic  basis  for  univarlant  isenthalpic  solidification  has  been  dis¬ 
cus.  along  with  evidence  for  its  occurrence  in  highly  supercooled  P4. 

2.  Univariant  Isenthalpic  solidification  appears  to  be  a  unary-system  analog  to  the 
solid-state  reactions  in  binary  alloys  known  as  "massive"  transformations. 

3.  Although  plane-front  solidification  at  constant  speed  is  theoretically  possible  in 
hypercooled  melts —  because  no  long-range  heat  transport  need  occur  — the  actual  mor¬ 
phology  of  the  solid/'llquid  interface  in  P4  remains  essentially  dendritic. 

4.  From  observations  on  a  macroscopic  scale,  it  appears  thai  a  well-defined  inter¬ 
face  sweeps  through  the  hypercooled  system  and  "massively"  transforms  the  melt  to 
solid.  However,  on  a  microscopic  scale  the  interface  probably  is  composed  of  protuber¬ 
ances  (scallops)  which  advance  by  release  of  latent  heat  Into  a  multi-dimensional  thermal 
diffusion  field  having  a  diffusion  length  on  the  order  of  only  10' 5  cm.  This  type  of  heat 
flow  is,  most  assuredly,  more  akin  to  dendritic  behavior  than  to  planar -front  behavior. 

5.  Solutions  to  the  time-independent  heat-flow  equation,  developed  by  Horvay  and 
Cahn  for  isothermal  dendrites,  we~e  used  in  the  present  analysis  with  some  modifica¬ 
tions  to  account  for  the  Gibbs-Thompson  effect  and  the  occurrence  of  appreciable  inter¬ 
face  supercooling.  Our  analysis  indicates  that  P4  solidities  in  modes  that  maximize  the 
Interface  supercooling  for  a  given  level  of  melt  supercooling.  Moreover,  the  data  indi¬ 
cate  that  uniform  attachment  kinetics  are  rate  controlling  when  P4  solidifies  with  a  nor¬ 
malized  interfacial  supercooling  greater  than  about  0.1.  The  linear  rate- constant  for  P4 
Is  estir  ated  to  be  17.7  ±  0.4  cm/sec  C°. 

6.  Layer-passage  limited  kinetics,  as  evidenced  by  the  appearance  of  faceted  growth 
morphologies  at  small  levels  of  supercooling  (Ad  <  1  C“),  change  over  to  uniform  attach¬ 
ment  kinetics  somewhere  in  the  range  of  melt  supercooling  between  1  C°  and  9  C°. 

7.  The  total  superc.  rHng  available  for  solidificatior  5s  distributed,  In  a  rather  com¬ 
plicated  way,  among  the  contributing  effects  of  molecular  attachment,  heat  flow,  and  in¬ 
terface  curvature.  No  one  effect  appears  to  be  extraordinarily  dominant  —  at  leas*  in 
terms  of  the  fractional  supercooling  expended  on  the  particular  effect  .wever,  at  large 
supercoolings  molecular  attachment  accounts  for  the  major  shai  e  of  tne  total  supercool¬ 
ing,  whereas  at  small  supercoolings  heat  flow  appears  to  be  the  largest  effect  as  long  as 
dendritic  growth  persists. 
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li.  MIM.CT 

Isenthalpic  solidification  of  a  pure  supercooled  liquid  is  shown  to  result  in 
either  a  two-phase  solid/liquid  mixture  in  invariant  equilibrium  or  a  single -phase, 
totally  solid  material  in  univariant  equilibrium,  depending  on  the  level  of  super¬ 
cooling  prior  to  solidification.  The  critical  supercooling  above  which  univariant 
equilibrium  is  obtained  is  large  for  metals  (hundreds  of  centigrade  degrees)  but 
much  smaller  for  certain  molecular  substances.  Experiments  on  white  phosphorus 
(alpha  P4)  show  that  the  critical  supercooling  (25.6  C°)  can  be  reached,  and  ex¬ 
ceeded,  easily.  Solidification  rate  measurements  taken  above  and  below  the  critical 
supercooling  for  P4  show  that  the  solid/liquid  interface  temperature  varies 
smoothly  with  melt  supercooling,  although  light-scattering  experiments  indicate 
that  rapid  changes  occur  in  the  extent  of  the  dendritic  zone  as  the  critical  super¬ 
cooling  is  approached  and  exceeded. 

A  method  for  extracting  interface  attachment  kinetics  from  solidification  rate 
data  was  examined  in  detail  and  applied  to  our  rate  measurements  on  P4.  We  find 
that  above  about  9  C°  supercooling,  P4  solidifies  with  linear  attachment  kinetics 
having  a  rate  constant  of  17.7  ±  0.4  cm/sec  C°.  Below  about  1  C°  supercooling,  P4 
solidifies  with  a  faceted  morphology  indicative  of  layer-passage  limited  kinetics. 
Between  1  C°  and  9  C°  supercooling,  transitional  growth  kinetics  occur.  These 
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results  are  in  qualitative  agreement  with  me  crystal  growth  theory  of  Cahn, 
et  al.,  which  predicts  that  attachment  kinetics  should  change  as  the  driving 
force  for  crystal  growth  is  varied  by  substantial  amounts. 


